*****************************************************************************
* Andrea L. Everett and Daniel C. Tirone					   				*
* Strategic Samaritanism Conditional Acceprtance at II		   				*
* Replication Do-File for Figures 2-6 in Manuscript and Appendices  		*
* For use with Everett-Tirone Replication Dataset II Conditional Acceptance	*
* STATA version 15															*
* To create figures hightlight from *START* to *END* and hit CTRL+D         *
*****************************************************************************


***Figure 2: Rope Ladder Plot, Aggregate Measures***
*START*
quietly xtreg f.lntotalaid3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, dydx(ae_che2 postche_1to3 postche_4to6) post
eststo total
quietly xtreg f.lntotalaid_prp10000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, dydx(ae_che2 postche_1to3 postche_4to6) post
eststo socgov
quietly xtreg f.lntotalaid_prp2000030000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, dydx(ae_che2 postche_1to3 postche_4to6) post
eststo economic
quietly xtreg f.lntotalaid_prp70000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, dydx(ae_che2 postche_1to3 postche_4to6) post
eststo humanitarian
coefplot (total, msymbol(D) label(ln Total Aid)) (socgov, msymbol(O) label(ln Soc & Gov Aid)) (economic, msymbol(T) label(ln Economic Aid)) (humanitarian, msymbol(S) label(ln Humanitarian Aid)), title(Marginal Effects and 95% CI) subtitle("Binary Emergency Indicators") xline(0) plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save Figure2, replace
*END*

***Figure 3: Predicted Aid, Millions of USD, Aggregate Measures***
*START*
quietly xtreg f.lntotalaid3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=0) post
quietly estimate store noche
quietly xtreg f.lntotalaid3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=1 postche_1to3=0 postche_4to6=0) post
quietly estimate store che
quietly xtreg f.lntotalaid3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=1 postche_4to6=0) post
quietly estimate store postche1to3
quietly xtreg f.lntotalaid3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=1) post
quietly estimate store postche4to6
coefplot (noche, msymbol(D) label(No Complex Emergency)) (che,  msymbol(O) label(Ongoing Complex Emergency)) (postche1to3, msymbol(T) label("1 to 3 years post-Complex Emergency")) (postche4to6, msymbol(S) label("4 to 6 years post-Complex Emergency")), vertical xlabel("") title("Total Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save F1, replace

quietly xtreg f.lntotalaid_prp10000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=0) post
quietly estimate store noche
quietly xtreg f.lntotalaid_prp10000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=1 postche_1to3=0 postche_4to6=0) post
quietly estimate store che
quietly xtreg f.lntotalaid_prp10000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=1 postche_4to6=0) post
quietly estimate store postche1to3
quietly xtreg f.lntotalaid_prp10000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=1) post
quietly estimate store postche4to6
coefplot (noche, msymbol(D) label(No Complex Emergency)) (che, msymbol(O) label(Ongoing Complex Emergency)) (postche1to3,  msymbol(T) label("1 to 3 years post-Complex Emergency")) (postche4to6, msymbol(S) label("4 to 6 years post-Complex Emergency")), vertical xlabel("") title("Soc & Gov Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save F2, replace

quietly xtreg f.lntotalaid_prp2000030000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=0) post
quietly estimate store noche
quietly xtreg f.lntotalaid_prp2000030000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=1 postche_1to3=0 postche_4to6=0) post
quietly estimate store che
quietly xtreg f.lntotalaid_prp2000030000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=1 postche_4to6=0) post
quietly estimate store postche1to3
quietly xtreg f.lntotalaid_prp2000030000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=1) post
quietly estimate store postche4to6
coefplot (noche, msymbol(D) label(No Complex Emergency)) (che, msymbol(O) label(Ongoing Complex Emergency)) (postche1to3, msymbol(T) label("1 to 3 years post-Complex Emergency")) (postche4to6, msymbol(S) label("4 to 6 years post-Complex Emergency")), vertical xlabel("") title("Economic Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save F3, replace

quietly xtreg f.lntotalaid_prp70000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=0) post
quietly estimate store noche
quietly xtreg f.lntotalaid_prp70000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=1 postche_1to3=0 postche_4to6=0) post
quietly estimate store che
quietly xtreg f.lntotalaid_prp70000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=1 postche_4to6=0) post
quietly estimate store postche1to3
quietly xtreg f.lntotalaid_prp70000ttl3avg ae_che2 postche_1to3 postche_4to6 lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 ae_che2=0 postche_1to3=0 postche_4to6=1) post
quietly estimate store postche4to6
coefplot (noche, msymbol(D) label(No Complex Emergency)) (che, msymbol(O) label(Ongoing Complex Emergency)) (postche1to3, msymbol(T) label("1 to 3 years post-Complex Emergency")) (postche4to6, msymbol(S) label("4 to 6 years post-Complex Emergency")), vertical xlabel("") title("Humanitarian Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save F4, replace

grc1leg F1.gph F2.gph F3.gph F4.gph, legendfrom(F1.gph) title("Binary Emergency Indicators") subtitle("Predicted Aid in Millions USD and 95% CI") plotregion(fcolor(white)) graphregion(fcolor(white)) iscale(.7)
graph save Figure3, replace
*END*


***Figure 4: Predicted Aid, Millions of USD, Violent Conflict Deaths***
*START*
quietly xtreg f.lntotalaid3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=(0(1)13) lnfdp_totalsource=0 lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Total Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono) 
graph save totalaid, replace

quietly xtreg f.lntotalaid_prp10000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=(0(1)13) lnfdp_totalsource=0 lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Soc & Gov Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save govcivaid, replace

quietly xtreg f.lntotalaid_prp2000030000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=(0(1)13) lnfdp_totalsource=0 lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Economic Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save econaid, replace

quietly xtreg f.lntotalaid_prp70000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=(0(1)13) lnfdp_totalsource=0 lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Humanitarian Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save humanitarianaid, replace

gr combine totalaid.gph govcivaid.gph econaid.gph humanitarianaid.gph, title("Violent Conflict Deaths") subtitle("Predicted Aid in Millions USD and 95% CI") plotregion(fcolor(white)) graphregion(fcolor(white)) iscale(.7)
graph save Figure4, replace
*END*

***Figure 5: Predicted Aid, Millions of USD, Displaced Persons***
*START*
quietly xtreg f.lntotalaid3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=(0(1)15) lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Total Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save totalaid, replace

quietly xtreg f.lntotalaid_prp10000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=(0(1)15) lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Soc & Gov Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save govcivaid, replace

quietly xtreg f.lntotalaid_prp2000030000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=(0(1)15) lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Economic Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save econaid, replace

quietly xtreg f.lntotalaid_prp70000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=(0(1)15) lnfdp_host=0) post
quietly marginsplot, ytitle("") title("Humanitarian Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save humanitarianaid, replace

gr combine totalaid.gph govcivaid.gph econaid.gph humanitarianaid.gph, title("Displaced Persons") subtitle("Predicted Aid in Millions USD and 95% CI") plotregion(fcolor(white)) graphregion(fcolor(white)) iscale(.7)
graph save Figure5, replace
*END*

***Figure 6: Predicted Aid, Millions of USD, Refugees Hosted***
*START*
quietly xtreg f.lntotalaid3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=0 lnfdp_host=(0(1)15)) post
quietly marginsplot, ytitle("") title("Total Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save totalaid, replace

quietly xtreg f.lntotalaid_prp10000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=0 lnfdp_host=(0(1)15)) post
quietly marginsplot, ytitle("") title("Soc & Gov Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save govcivaid, replace

quietly xtreg f.lntotalaid_prp2000030000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=0 lnfdp_host=(0(1)15)) post
quietly marginsplot, ytitle("") title("Economic Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save econaid, replace

quietly xtreg f.lntotalaid_prp70000ttl3avg lnucdp_allconflictdeaths lnfdp_totalsource lnfdp_host lnemdat_total_affected lnGDPpc_constant2005USD lntrade_pctgdp lninfant_mortality lnpopulation_total FH_Combined PctAgreeUS year_count if lnfdp_totalsource!=.&year>=1989&oecd_highincome!=1, fe vce(cluster recipient_cowcode)
quietly margins, expression((exp(predict(xb)))/1000000) atmeans at(lnemdat_total_affected=0 lnucdp_allconflictdeaths=0 lnfdp_totalsource=0 lnfdp_host=(0(1)15)) post
quietly marginsplot, ytitle("") title("Humanitarian Aid") plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono)
graph save humanitarianaid, replace

gr combine totalaid.gph govcivaid.gph econaid.gph humanitarianaid.gph, title("Refugees Hosted") subtitle("Predicted Aid in Millions USD and 95% CI")plotregion(fcolor(white)) graphregion(fcolor(white)) scheme(s2mono) iscale(.7)
graph save Figure6, replace
*END*
